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Introduction - Digital Filters (INCT) 


音声 や 画像 な ど を パソ コン な どの 計算 機 上 で 扱う 場合 、 サ ンプ リ 
ング され た 信号 と し て 取り 扱う こと に な る と の よう な 形態 で 信号 
を 取り 扱う 手法 は ディジ タル 信号 処理 技術 と し て 知ら れ て いる 
現在 で は ヽ 情報 通信 の 分 野 は も ちろ ん ヽ 広 く 工 学 分 野 に お いて 、 ヽ 
ディ ジタル 信号 処理 に 関す る 理解 が 欠か せな く な っ て き て いる 


そこ で ヽ ゝ 本 テー マ で は ディ ジタル 偽 号 処理 技術 の 基礎 の 一 つ で あ 
る ディ ジタル フィ ル タ を 理解 する こと を 目的 と する 具体 的 に は 、 
イン パル ス 応 答 か ら デ ィ ジ タル フィ ル タ の 特性 を 計算 する 方 法 や 、 
信号 の フィ ル タ リ ング 処理 ( 辿 み 込み 演算 ) 手 順 を 理解 し 、 プ ログ ラ 
ム の 実装 を 行う 


Overview ot Digital Filters (INCT) 


線形 シフ ト 不 変 シ ステ ム と イン パル ス 応 答 


ディ ジタル 信号 処理 を 実行 する 処理 プロ ッ ク の モデ ル と し て ヽ 図 1 
に 示す ブロ ッ ク 図 を 考え る < 


入力 信号 xn) | シス テム | 出力 信号 y(n) 


シス テム の 入出 力 関係 


図れ 1 で ヽ 、 z(z) お よび ヶ ( ヵ ) を それ ぞ れ 入力 倍 号 お よび 出力 信号 と す 
る 。 ここ と で ヽ g(%) か ら ?(z) を 生じ る 処理 プロ ッ ク を ヽ シ ステ ム と 
呼 信 - と に する * シス テム が 線形 性 、 シ フト 不変 性 [1] な どの 条件 を 
満た す と き ヽ と の シス テム を 緑 形 シ フト 不変 シス テム と いう ま 
た ヽ シ ステ ム を 実 時 間 処 理 で 実現 する に は 、 ヽ 因果 性 [1] が 要求 され 
る 本 テー マ で 扱う シス テム は すべ て 因果 性 を 満た す も の と する 


線形 シフ ト 不 変 シ ステ ム に ディ ジタル 単位 イン パル ス 6(? ヵ ) 
Edquation: 


を 入力 し た と き の 出 力 ヵ ( ヵ ) を 、 イ ン パ ルス 応答 と 定義 むる > イン 
パル ス 応 答 は ヽ ュ シ ステ ム の 性 質 を 示す 重要 な 値 で ある 


由 み 込み 演算 と ディ ジタル フィ ル タ 


ある シス テム が 線 形 性 と シフ ト 不 変性 を 満た す 場 合 、 そ の シス テム 
の イン パル ス 応 答 ヵ ( ヵ ) が 選 え られ れ ば ヽ 任意 の 入力 z( ヵ ) に 対す る 
出力 y(z) は ヽ 式 (2②) に より 計算 する こと が で きる 


Equation: 


7) = ニ に -oo (6 のみ 一 
ー クン に -oo (が)z( み 一 


さら に ヽ ゝ ある 入力 信号 に 対し て 晴 み 込み 演算 を 実行 する こと に よ 
り ヽ 入力 信号 か ら あ る 特定 の 周波 数 成分 だ け を 抜き 出す な どの フィ 
ル タ 効 果 が ある と と が 知 られ て いる 加え て この フィ ル タ 効 果 は 、 
イン パル ス 応 答 を 変化 させ る こと と に より 様々 に 変化 させ る こと が で 
きる この こと と か らい ヽ 晴 み 込み 演算 を 実行 する シス テム は ヽ デ ィ ジ 
タル フィ ル タ と も よ ば れ て いる イン パル ス 応 答 が 有限 時 間 だ け 持 
続 す る シス テム を FIR(Finite Impulse Response) シ ステ ム ( ま た は FIR 
フィ ル タ ) ヽ 無限 時 間 持 続 す る フィ ル タ の こと を IIR(Infinite Impulse 
Response) シス テム (また は IIR フ ィ ル タ ) と 分 類 し て いる ・ 本 テー マ で 
は ヽ EIR フ ィ ル タ を 扱う こと と する 


イン パル ス 応 答 と 周波 数 特性 


以上 の よう な フィ ル タ 効 果 を 弓 を える シス テム の 特性 を 、 そ の シス 
テム の 周波 数 特性 と し て 表現 する こと が で きる 。・ イ ン パ ルス 応 管 
ヵ ( ヵ ) を 有する シス テム の 周波 数 特性 万 (e) は ヽ ゅ を 正規 化 角 周 波 
数 と し て ヽ 次 式 で 定義 され る 
Equation: 


万 (e7) = ゞ 、 (7)e 7 


77ーOO 


シス テム が EIR フ ィ ル タ の 場合 、 イ ン パ ルス 応答 ヵ ( ヵ ) は 有限 時 間 
し か 持続 し な い ぃ またい ヽ 因 果 性 を 満た す サシ ステ ム の 場合 因果 性 の 
定義 より ヽ ヵ ( ヶ ) = 0( ヵ < 0) で ある の で ヽ 因 果 性 を 満足 し イン パル 
ス 応答 長 お が WW で ある FIR シ ステ ム の 周波 数 特性 は 次 式 と な る * 


Edquation: 


MV の と と を フィ ル タ 長 と いう こと も ある 


周波 数 特性 の 計算 


周波 数 特性 に は ヽ 、 そ の 大 き さ を 表す 振幅 特性 と ヽ 入 出力 信号 間 の 
位相 の 変化 を 表す 位相 特性 と が ある 。 ま た ヽ 位 相 特 性 を 別 の 形 で 表 
現し た 特性 と し て ヽ 群 遅延 特性 が ある < 


④ 式 より ヽ 振 幅 特性 互 (e2) は ヽ 以下 の よう に 実 部 と 虚 部 に 分 け て 
表す こと が で きる ・ 


Edquation: 
-1 -1 
万 (e) = ん (7) cos o7 十 7 - の (7) sin a| 
7 王 0 7 三 1 


こと で ヽ 右 辺 第 一 項 は 実 部 右辺 第 二 項 は 虚 部 で ある ・ よ っ て 振幅 
特性 互 (e7?) の 絶対 代 は ヽ (⑤5) 式 の 絶対 債 を と っ て 以下 の よう に 計算 
0 
Equation: 


ge の )| = V+ 


式 (6) で 刀 。 お よび 了 万, は 、 ヽ それぞれ (5) 式 に お ける 実 部 と 虚 部 を 表す 
も の と する > 振幅 特性 の 大 き さ は ヽ し ば し ば [dB] で 表示 され る こ 
の 場合 の 振幅 (ampitude) [dB] は ヽ 常用 対数 を と っ て 次 式 で 表 さ れる < 
Edquation: 


20 1og。 | 万 (e が の ) | 


⑦ 式 で ヽ 正規 化 角 周波 数 。 を 細か く 変 化 さ せ て を 横 軸 ヽ 振幅 特性 
を 縦 還 と し て グラ フ に 示せ ば ヽ ゝ 振幅 特性 の グラフ を 得る こと が で き 
る 。 ま たい ヽ 骨 波数 を 正規 化 周波 数 ナー= 愛 で 表現 する こと も 多い < 
正規 化 周 小数 げ は ヽ 周波 数 刀 [Hz]j を サン プリ ング 周波 数 戸 [Hz] で 
割っ た 値 で ある 

位相 特性 9(。o) は ヽ 振幅 特性 万 (e め ) を 極座標 表示 し た と き の 備 角 で 
ある < 数 値 計算 に より 9(。) の 値 を 求め る 場合 に は ヽ 次 式 を 用 いる 
こと が で きる (た だ し ヽ そ の 処理 系 で の 逆 正 接 関数 一 た と えば atan( 
人 の 主 値 の 範囲 に つい て 理解 し て お く こ と が 必要 で あろ 

ブ O 


Equation: 


9() =tan「~ | 


と 


さら に ヽ 群 遅延 特性 7 の 定義 は ヽ 位相 特 性 の 傾き と し て 次 式 で 避 え 
られ る 


Edquation: 


[9o7zp/e| 


以上 を まとめ る と ヽ イ ン パ ルス 応答 ヵ (z) の 値 が 既知 な ら ば ヽ 周波 数 
特性 (振幅 、 位相 ヽ 群 遅延 特性 ) を 計算 で きる こと に な る < ある イン 
パル ス 応 答 に つい て ヽ 実 際 に 周波 数 特性 を 計算 し た 例 を 図 2 に 示 

す 図 2 ぬ ) が イン パル ス 応 答 で ヽ 図 20b) か ら (① が それ ぞ れ われ 振幅 、 位 
相 ヽ 群 遅延 特性 で ある < 
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周波 数 特性 の 計算 例 


伝達 関数 に つい て の 補足 


式 (2) を z 変 換 す る と ヽ 次 式 と な る < 
Equation: 


上 式 の 互 (>) を 伝達 関数 と 呼ぶ * 伝達 関数 万 (>) は インパ ルス 応答 
ヵ () の ヶ 変 換 で あり ヽ 、 デ ィ ジ タル フィ ル タ の 入出 力 関係 を z 領域 


で 表し た 式 で ある ゝ 式 (3) と 式 (11) を 比較 する と ヽ 晴 み 込み 演算 が > 
領域 で は 積 で 表現 で きる こと が わか る 


な AS 上 E 民 に ッ ー ン の を 作る らら 次 式 ど な る ? 
Edquation: 


Y(e の ) 互 (e%) 双 (e の ) 


え (e?) お よび 予 (e?) は 、 ヽ それ ぞ わ れ 人 入力 信号 z(z) お よび 出力 信号 
9(?) の 離散 時 間 フ ー リ エ 変 換 (周波 数 スペ ベクトル ) で ある ・ 式 (11) は ヽ 
畳み 込み 演算 を 周波 数 領域 で 表現 し た 式 に 相当 する < 


Analysis and Implementation of Digital Filters - (INCT) 


[課題 1] 下図 に 示す イン パル ス 応 答 を 有する ヽ イ ン パ ルス 応 
答 長 の EIR シ ステ ム を 考え る 
0.6 


0 5 10 15 20 
Time in Sample 


イン パル ス 応 答 (N=3) 


この と きい ヽ 以 下 の 手 順 で この シス テム を 解析 し ヽ 振 幅 特 性 を 計算 せ 
よ > また ヽ 上 図 の 1 インパ ルス 応答 を 持つ FIR デ ィ ジ タル フィ ル タ と 
し て 実現 せよ 


(1) シ ステ ム の ブロ ッ ク 図 を 示せ 


(2) イ ン パ ルス 応答 か ら シ ステ ム の 振幅 特性 を 計算 し 、 グ ラフ 
と し て 表示 せよ ・ た だ し ヽ 横 軸 は 正規 化 周波 数 ヽ 縦 還 は [dB] と す 
る o 


(3) 信 号 デ ー タ (ファ イル 名 : [missing_resource: signal.txt] どじ て 
える は ヽ 100 個 の サン プル 値 か ら な る 信号 で ある (Figure 2) * 
計算 機 上 の グラ フッ ー ル を 用 いて この 依 号 を 図示 し た と き ヽ Figure 2 
の よう に な る こと を 確認 せよ (Figure 2 は 、 サ ンプ ル 値 を 強調 する た 
め ヽ 本 グ ラフ 状 に 表示 し て いる が ヽ 折線 グラ フ で 表示 し て も よ 

い )・ 


(4)x(n) を 入力 信号 と し て この シス テム に 加え ヽ EIR フ ィ ル タリ ング を 
行なっ た 出力 ym) を 求め る プロ グラ ム を 作成 せよ 


⑤5) 6 2 の よう な 形式 で 図示 し ヽ 入力 信号 と の 波形 の 変化 
を 考察 せ 
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入力 と し て 使用 する サン プ 


ル 値 信号 
以上 を 含む 内 容 を ヽ イ ン パ ルス 応答 を 


有する FIR フ ィ ル タ の シミ ュ レ ーション 結果 と し て 整理 し て 報告 せ 
中 o 


た だ し ヽ ゝ 結果 の 報告 の 方 法 に は 注意 し ヽ 自 ら が 問題 意識 を 持っ て 
実験 を 行ない 報告 する ヽ と いう 立場 か ら 内 容 を 整理 し て 報告 する こ 
と 5 


・ 何 を 目的 と し て ヽ ど の よう な 内 容 の 実験 を 行 な お うと し て いる 
レポ ボート 


・ ど の よう な ブロッ ク 了 図 の シス テム で ヽ ど の よう な 処理 を 実行 し 
て いる の 


・ ど の よう な 信号 を 加え ヽ 、 ど の よう な 処理 結果 を 得 た の か ゃ 
の 結果 ヽ 、 ど の よう な こと が 育 え る の か ゃ * 


上 記 の よう な 点 に 注意 を 払い ヽ 報告 書 を 作成 の こと < 


[課題 2] 以下 の 図 に 示す インパ ルス 応答 (ファ イル 名 : 
[missing_resource: h10.txtl 
0.6 
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イン パル ス 応 答 (N=10) 


[課題 3] 課題 2 と 同様 に 、 以下 の 図 に 示す イン パル ス 応 答 ( ファ イル 


名 : [missing_resource: h21.txt] 
06 
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イン パル ス 応 答 (N=21) 


音声 斉 号 (男声 ) に 周波 数 4[kHz] の 正弦 流 音 号 が 妨害 斉 号 と し 
れ た Figure 5 に 示す 信号 が ある この 信号 か ら 妨 害 依 号 を 取 
音声 信号 の み を 取り 出し た い * 10 ィ ル タ を 設計 
まで の 課題 で 作成 し た た た み 込 み 計算 プロ グラ ム ( フ ィ ル 


人 


タリ ング プロ グラ ム ) に より 妨害 信号 を 取り 除い た 信号 を 生成 せよ 
な お ヽ こ の 信号 の サン プリ ング 周波 数 は Fs=11.025[kHz] で ある < 

※ な お ヽ こ の 課題 に つい て は 追加 資料 を 配布 し ます の で ヽ ゝ 参考 に し 
を ド さ ツル ュ o 
20000 


20000 40000 
サン プル 時 刻 


正弦 波 信号 が 妨害 信号 と し て 付加 され た 音声 信 
号 流 形 (Fs=11.025[kHz]) 


Assignments - Digital Filters (INCT) 


(① 図 4 の 入力 信号 は ヽ 正弦 渡 に 雑音 を 加え た 依 号 を サン プリ ング し 
た 信号 で ある 図 4 の 横 軸 を よく みる と わか る よう に ヽ こ の 信号 は 正 
弦 波 1 周期 を 20 点 で サン プリ ング し て いる この 正弦 波 の ディ ジタル 
周波 数 (jE 規 化 周波 数 ) を 示せ < 


(2) 各 課題 で の FIR フ ィ ル タ の 振幅 特 隆 か ら ヽ EIR フ ィ ル タリ ング 前 後 
の 信号 波形 の 変化 を で きる だ け 定 量 的 に 説明 せよ (x(m) と y(m) に つい 
て ` ヽ 信号 の 天 き さき 。 波形 の 変化 を 比較 する ) 


(3)FIR フ ィ ル タ と IIR フ ィ ル タ の 特徴 を 調べ ヽ 比 較 せ よ * 


(4) イ ン パ ルス 応答 の 値 を さま ざま に 変え る こと で い ヽ 特 性 の 異な る 
フィ ル タ を 実現 で きる 希 思 の 処理 を 行う た め に 適切 な イン パル ス 
応答 の 値 を 決定 する こ と が ヽ デ ィ ジ タル フィ ル タ の 設計 に 相当 す 
る < どの よう な ディ ジタル フィ ル タ 設 計 法 が あり ヽ 、 そ れ ぞ れ ど の よ 
うな 特徴 を 持つ か を 調べ て みよ 


5) 感想 を 記せ 。 


AdVanced Problems - Digital Filters (INCT) 


(1 イン パル ス 応 答 か ら 位 相 特 性 を 計算 する プロ グラ ム を 作成 せよ * 
た だ し ヽ tan~~ 全 が 数 学 的 に 未定 義 で も ある と と に 注意 し て プロ グ 
ラム を 作成 し な いと ヽ 振 幅 の 値 が せ ロ に 近づく 周波 数 で アプ ログラム 
が 実行 時 エラ ー を 起こ す 場 合 が ある * 


(2) イ ン パ ルス 応答 か ら 群 遅延 特性 を 計算 する プロ グラ ム を 作成 せ 
よ ・ こ の 結果 か ら ヽ ( ) と ( ) の 波形 の 時 間 的 な ずれ を 定量 的 に 
考察 せよ 


(3) 位相 特性 と 群 遅 延 特 性 は ヽ 実質 的 に 同じ こと を 異な る 表現 で 示し 
て いる に すぎ な いこ と を い ヽ 実 験 結果 か ら 考 察 せ よ 


(4)http:/momijii.ishikawa-nct.ac.jp に より ヽ ゝ 各種 の FIR お よび IIR デ ィ ジ 
タル フィ ル タ が 設計 可能 で ある 様々 な フィ ル タ の 設計 を 行なっ て 
みよ また ヽ 得 られ た フィ ル タ 係 数 を 利用 し て 、 ヽ ディ ジタル フィ ル 
タ を 実現 せよ < 


